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Abstract 

In this letter, we show that the formulae of Bray and Moore for the average 
logarithm of the number of metastable states in spin glasses can be obtained by 
calculating the partition function with m coupled replicas with the symmetry 
among these explicitly broken according to a generalization of the 'two-group' 
ansatz. This equivalence allows us to find solutions of the bm equations where 
the lower 'band-edge' free energy equals the standard static free energy. We 
present these results for the Sherrington-Kirkpatrick model, but we expect 
them to apply to all mean- field spin glasses. 
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There has been recently a renewed interest in the study of metastable states in spin 
glasses. Current problems range form their potential relevance in spin glass dynam- 
ics to the formulation of renormalization group analysis in the presence of many 
minimas of the action. Nevertheless, it is interesting to use our present knowledge 
to shed some light old problems and to understand the connection between new and 
older techniques. 

One recent development is the use of coupled real replicas as a way to explicitly 
break ergodicity revealing the 'complexity' of the system Jl], It was not known 
up to now whether or not these techniques were equivalent to the now classical 
computations of the number of tap solutions of Bray and Moore (bm) || |]]. It 
has been shown to be true for the spherical p-spin model, the Potts glass and for 
part of the spectrum of the random orthogonal model (ROM) || ^ [/J. A much 
older question that remained partly unanswered is whether the computation of BM 
for the Sherrington-Kirkpatrick model (sk) using continuous RSB gives at the lower 
'band-edge' the same free energy as the static calculation. It was shown in Ref. 
that they agreed near the transition, at least to the first non-trivial order in the 
reduced temperature t. The BM computation involves a function Q(x) similar to 
the usual Parisi function q(x) but in this computation the critical value of x (x c ) 
is equal to ~ near T c while its static counterpart only equals 2t. Finally, a never 
explained coincidence is the fact that the 'direct average' of BM both for SK and ROM 
can be obtained by calculating the partition function under the two-group replica 
symmetry-breaking scheme. If one does not let the number of replicas n go to zero, 
the two-group partition function gives the Laplace transform of the number of TAP 
solution with n conjugated to —(3f |7|, H]. 

The contribution of this letter is to show that a generalization of the two-group 
ansatz makes the connection between the 'real replicas' and the BM framework. This 
connection in turns helps us understand why, in the 'naive' two-group approach, n 
becomes the Laplace conjugate of the free energy. The two-group model, in its 
unbroken form, gives us the proper ansatz to look for solutions of the BM saddle 
point that reduces to the static free energy at the lower band-edge. All these points 
are illustrated on the SK model but we believe that they should apply to any mean- 
field spin glass. Unfortunately, a crucial point stills eludes us, that is the physical 
significance of the two-group ansatz and the difference, if any, between metastable 
states that break the two-group symmetry and those that don't. 

In Ref. [U it was argued that the number of pure states in a glassy system can be 
computed by introducing m real replicas of the system. More precisely, for a mean- 
field model with quenched randomness where one introduces n replicas to average 
over the disorder, we have the following relation 

max [(log N{f)) - m(3f] = lim - log (Z mn ) , (1) 

/ n.->0 n 

where Af(f) is the number of metastable states as a function of their free energy 
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and Z mn is the partition function ofnxm copies of the original system with the 
replica symmetry explicitly broken into n groups of m replicas. 

It is useful to consider both / and (logA/"(/)) as function of the independent 
parameter m through the relation / = f*(m) i.e. the particular / that maximizes 
the l.h.s. of Eq. ([!]). Varying m then traces out the (log jV(/)) vs / curve, in general 
a bell shaped curve (e.g. SK or ROM || or an 'interrupted bell' shape (e.g. p-spin, 
both spherical and Ising [||, ||). The maximum of this curve is given by m = 
(white average) while increasing m favors states of lower free energy. In particular 
m = 1 picks out the states that contribute to the Gibbs measure. Finally, we call 
the lower band-edge the infimum of the free energies for which (logA/"(/)) > 0. 

This formula along with the simplest replica symmetry ansatz gives a very good 
description of the metastable states for mean-field models with 'first-order rsb' 
(e.g. spherical p-spin, the Potts glass and ROM). It predicts that at the dynamical 
transition, a large number of metastable states start to contribute to the free energy 
(extensive complexity) and that at the RSB transition the complexity ceases to be 
extensive [|IJ. In the case of SK, a simple one-step ansatz in Eq. ([]]) always gives 
either zero or a negative result for the complexity. As we shall see, the presence 
of metastable states in SK can be uncovered by a more complex replica symmetry 
breaking scheme. 



The partition function of the replicated SK model is given by (see e.g. [10 

P 2 



(Z mn ) = Ext exp { N 
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mn - TrQ 2 ) + log £ exp ±/? 2 £ Q>X 



, (2) 



where the extremum is taken over matrices Q^ b (1 < a,b < n and 1 < i,j ' < m) 
with zeros on the diagonal Q™ a = 0. 

The point of this letter is to show that the formula for the average logarithm of 
the number of tap solution of Ref. |4[] can be obtained from Eq. (P using a direct 
generalization of the 'two-group' replica symmetry breaking |§. We will use the 
following ansatz on each of the n 2 sub-matrices Q a f, of size m x m (i.e. the matrices 
Q % l h with (a, b) fixed and (i, j) varying form 1 to m): 

V rn-y 

where y is to be taken to infinity. A convenient way to parameterize the matrices 
Qab and Q~ b in order to obtain a finite result as y — ► oo is 

<& = Qo6±^ + %. (4) 

y y 

Note that the 'diagonal matrices' Q aa have the extra constraint that their diagonal 
elements Q™ a be zero; their off-diagonal elements, in general non-zero, are also pa- 
rameterized according to Eq. (M. We do not need to specify the structure of the 
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three matrices Q ab , A ab and B ab in order to make the connection with the BM for- 
mulation, but in practice they will be taken to have the usual utrametric structure 
(see e.g. IfLOt ). Using this ansatz, we can compute the partition function (H) in a 
manner similar to that of Ref . pi] . We find 



(Z mn ) = Extexp {nN [m\og2 + \f3 2 m{l - Q) 2 - (3 2 A{1-Q) 



-N(3 2 



\Kb + QabBab + mQ ab A ab + \m 2 Q 2 ab 

■ ab 



+ JVlog/L (5) 



where 



1 = /In^£lI(T?^)-p|-i:[ Io tanh- mo 

+ \m log(l - m 2 a ) +P 2 \QabX a Xb + B ab m a m b + A ab m a x b \ , (6) 

ab J 

and where A and Q are the diagonal terms A aa and Q aa , note that the sum over 
(a, b) also includes diagonal terms. If we put this formula in Eq. ([]]), we find that 
the average log number of metastable states is identical to the average log number 
of tap solutions as computed in Ref. [f|. We recall that in that paper, the number 
of tap solutions is written as the integral over all configurations of an appropriate 
delta function, replicas are introduced to average the logarithm over the disorder. 
The equivalence between our formula and Eq. (11) of Ref. ||] can be rendered more 
obvious by the following change of variables, 



dm„ 



(5 2 (2B ab + mA ab + \m 2 Q 



Vab = (3 2 Qab , Pab = (A ab + \wQ 

q = Q, A = (3 2 (A+ \mQ) , A = (3 2 (B + \mA + \m 2 Q) and 



ab 



u 



-m . 



(7) 



and a shift of the integration variable 

In Ref. § it was shown that the original two-group partition function (with n 
in general non-zero) gives the Laplace transform of the annealed averaged number 
of metastable states with n conjugated to —ftf. We recover this annealed result 
by setting the 'off-diagonal' matrices Q ab to zero, leaving us only with n diagonal 
blocks ofraxm two-group matrices. This explains the strange scaling with n, i.e. 
the size of the total matrix Q^ b has to be taken to zero, while the ij block size m 
remains finite and becomes the Laplace conjugate of the free energy. Similarly, using 
an ij-two-group over an RS Q ab , we recover the 'innocent replica' treatment, in this 
case, the full matrix Q^ b can be viewed as a one-step matrix with blocks of size m 
within which lies a two-group structure. 

Our formulation of the BM problem allows us to find a class of simple solutions 
to the saddle point equations, i.e. the 'unbroken two-group'. Indeed one can show 
that A ab = B ab = forms self-consistent solution of these equations. The problem 
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is then reduced to finding the extremum with respect to a single matrix Qa6 with a 
diagonal term Q. The computation with a fc-step RSB ansatz on the matrix Q a ^ is 
equivalent to the usual computation of the partition function with k + 1-steps with 
RSB parameters < mx\ < . . . < mxt < m < 1. In particular, since Eq. (|I|) implies 

log (A/7/)) = -m 2 ^^ when (3f = J-( m F(m)) , (8) 
dm dm 

where 

F(m) = lim — log (Z mn ) , (9) 
n-*o 

the condition for the lower band-edge is just the usual minimization with respect to 
the last RSB parameter m. 

In the simplest case, the 'innocent replicas' admit a solution that is equivalent 
to the one-step partition function. We have explicitly checked that the expansion 
of (pi) near T c (Eqs. (12,13)Q of Ref. §) admits a saddle point of the form (|7|) with 
A = B = 0, we find, for the lower band-edge 

q = t + § t 2 + 0(t 3 ) , q D = §t + 0(t 2 ) , m = ft + 0(t 2 ), 

and 

Pf/N = - log 2 - \{? + |t 3 + iit 4 + (5099/7290)t 5 + 0(t 6 ) . (10) 

This solution should be compared with that of Ref. ||, where it was found q = 
0.3471t + 0(t 2 ). Our solution can be followed to other values of the parameter m, 
but we find it to be only valid at the lower band-edge. As expected, for m > |t, it 
gives a negative log-number of tap solutions. Whereas for m < |t, / decreases with 
decreasing m indicating that we are on an unphysical branch. The correct solution 
for m < |t must therefore include two-group symmetry breaking. 

The same 'unbroken two-group solution' can be found for continuous RSB near 
T c . We find a solution very similar to that of Ref. 0], except for the value of m 
which in turns changes the value of all other parameters. We find 

q = t + t 2 + 0(t 3 ) , q(x) = tx(x<x c = l), m = 2t + 0(t 2 ) , 

and 

(3f/N = - log 2 - \(5 2 + |t 3 + iit 4 + l t> + 0(t 6 ) (11) 

The apparent discrepancy between the value of x c — 1 and that of the usual static 
solution (x c = 2t) is understood by realizing that it is not the matrix Q a b which 
should equal the static one but rather the full matrix Q^ b . When n goes to 0, the 
full matrix can be parameterized by a function q{ n n(x), 

qM\(x) = q{x/m) , < x < m\ q {u \\(x) = Q , m < x < 1 , (12) 

in perfect agreement with the static solution. 

1 We think that the factor (1 — 4Q) in that second equation should read (1 — 5Q). 



5 



Acknowledgments 

We wish to thank M. Ferrero and R. Monasson for useful discussions. 



References 

R. Monasson, The structural glass transition and the entropy of the metastable 
states, preprint |cond-mat/9503166| . 

S. Franz and G. Parisi, Recipes for metastable states in spin glasses, preprint 
|cond-mat/9503167| . 

A. J. Bray and M. A. Moore, Metastable states in spin glasses, J. Phys. C: 
Solid St. Phys. 13 (1980) L469. 

A. J. Bray and M. A. Moore, Metastable states in the solvable spin glass model, 
J. Phys. A: Math. Gen. 14 (1980) L377. 

A. Crisanti and H.-J. Sommers, On The TAP Approach to the Spherical p-spin 
SG Model, preprint TNT-94-5, |cond-mat/9406051| . 

T. R. Kirkpatrick and P. G. Wolynes, Stable and Metastable States in Mean- 
Field Potts and Structural Glasses, Phys. Rev. B. 36 (1987) 8552. 

G. Parisi and M. Potters, Mean- field equations for spin models with orthogonal 
interaction matrices, preprint |cond-mat/9503009 . 



A. J. Bray and M. A. Moore, Broken replica symmetry and metastable states 
in spin glasses, J. Phys. C: Solid St. Phys. 13 (1980) L907. 

H. Rieger, The Number of Solutions of the Thouless-Anderson-Palmer Equa- 
tions for the p-Spin- Interaction Spin Glasses, Phys. Rev. B 46 (1992) 14655. 

M. Mezard, G. Parisi, M. A. Virasoro, Spin glass theory and beyond, World 
Scientific (Singapore) 1987. 



6 



